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Abstract 

Athermal lattice gases of particles with first neighbor exclusion have been studied for a long 
time as simple models exhibiting a fluid-solid transition. At low concentration the particles occupy 
randomly both sublattices, but as the concentration is increased one of the sublattices is occupied 
preferentially. Here we study a mixed lattice gas with excluded volume interactions only in the 
grand-canonical formalism with two kinds of particles: small ones, which occupy a single lattice 
site and large ones, which occupy one site and its first neighbors. We solve the model on a Bethe 
lattice of arbitrary coordination number q. In the parameter space defined by the activities of both 
particles. At low values of the activity of small particles (zi) we find a continuous transition from 
the fluid to the solid phase as the activity of large particles (^2) is increased. At higher values of 
zi the transition becomes discontinuous, both regimes are separated by a tricritical point. The 
critical line has a negative slope at zi = and displays a minimum before reaching the tricritical 
point, so that a reentrant behavior is observed for constant values of Z2 in the region of low density 
of small particles. The isobaric curves of the total density of particles as a function of zi (or Z2) 
show a minimum in the fluid phase. 
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I. INTRODUCTION 



In the beginning of statistical mechanics, the interest was focused mainly on fluids, and 
the pioneering work on phase transitions concentrated on the liquid-gas transition, an exam- 



ple of this kind is the van der Waals equation of state In his phenomenological theory, 
van der Waals considered the effect of hard core excluded volume interactions and attractive 
interactions in the thermodynamic behavior of the fluid. While attractive interactions be- 
tween molecules are essential to produce the liquid-gas transition, it was later reahzed that 
even if only hard core interactions are considered, interesting effects arise in the models. 
Since in this case all allowed microscopic configurations of the system have the same energy, 
such models are athermal, and all thermodynamic effects are of entropic origin. Much is 
known about the continuous version of these models, usually called hard sphere systems. 
They were studied by a variety of techniques and a fluid-solid phase transition is found. 
t is worth recalling that in the pioneering work on the Monte Carlo simulational procedure 
, the physical system studied was a fluid of hard disks. 

The fluid model with excluded volume interactions only may also be defined on a lattice, 
so that the positions occupied by the particles are restricted to sites of a lattice. In the 
simplest version of such models, the only constraint is that if a site is occupied by one 
particle, no other particles may be placed on it. In the grand-canonical ensemble, this model 
reduces to the Ising model without the interaction term, and is therefore trivially solved. 
No singularities are found in the thermodynamic functions, as expected. More interesting 
results are obtained if the range of excluded volume interactions is increased. This leads to 



a variety of models, and we refer to a recent paper where some o 



these models were studied 



a. 



If a particle placed on 



using simulations for a comprehensive survey of the literature 
a site excludes this site and its first neighbors, indeed a fluid-solid continuous transition is 
found. For bipartite lattices such as the hypercubic lattices, at lower densities, the particles 
occupy the lattice sites randomly, but above a critical density one of the two sublattices is 
occupied preferentially by the particles, so that the order parameter may be defined as: 

ij = 2{\pA-pB\), (1) 

where pa,b is the number of particles in sublattice A or B, respectively, divided by the 
number of sites in the lattice, thus assuming a maximum value equal to 1/2. This model 
has a long history, it has been mentioned in the classical review by Domb (5| , and has been 



studied using the virial expansion and the Bethe approximation by Burley shortly after j6|. 
Since then, the thermodynamic behavior of the model has been investigated using a variety 
of analytical and numerical methods, which indicate a continuous phase transition in the 
Ising universality class. On the square lattice, precise estimates for the critical chemical 
potential ~ 1.33401510027774(1) and the density of particles at the transition ~ 
0.3677429990410(3) were obtained using transfer matrix and finite size scaling extrapolation 
techniques, and the Ising edtieal exponents were verified with high precision Q A recent 
simulational study of this model on the square and cubic lattices may be found in [Sj]. On the 
triangular lattice, this model is known as the hard hexagon model, and was solved exactly 
by Baxter js]. A continuous transition was found at the critical activity Zc = exp(/ic) = 
(11 + 5-\/5)/5. The critical exponents are in the 3-state Potts universality class, as would be 
expected considering that in the high density phase the hexagons occupy preferentially one 
of the three sublattices. 

An interesting generalization of the model is to consider a gas with both small and large 
particles, where the small particles occupy a single site and the large ones the site and its 
first neighbors. Let the activity of small (large) particles be zi {z2). On a square lattice, we 
may represent both particles as squares, as is shown in Fig. [H Since a continuous transition 
occurs when only large particles are present and no transition is found for small particles 
only, one may ask what happens in intermediate situations, with both types particles on the 
lattice. This model was studied on the square lattice by Poland HI] using high density series 
expansions in the grand-canonical formalism, and evidence was found that a tricritical point 
should be present in the phase diagram. In this paper we solve the model on a Bethe lattice of 
arbitrary coordination number g > 2 in the grand-canonical ensemble. As the fugacity of the 
small particles is increased starting from zero, the fluid-solid transition remains continuous 
up to a certain value, above which it becomes discontinuous. Therefore, we find a tricritical 
point in the phase diagram of the model. 

It is worth mentioning that a slight modification of the two-particle model on the square 
lattice makes it exactly solvable in a particular case, by allowing it to be mapped on the 
Ising model for which the exact solution is known for zero magnetic field. This particular 
case corresponds to {1 + Zi^ = Z2 [3]. In this lattice gas, the large particles are placed 
on the centers of the elementary squares of the lattice defined by the sites where the small 
particles are located, as is illustrated in figure El Both models are equivalent in the absence 
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of small particles, but we notice that when = we know the exact solution only for 

Z2 = l. 

In section [IT] we define the model and solve it on the Bethe lattice. The thermodynamic 
properties of the model are discussed in section IIIH and final discussions and the conclusion 
may be found in section IIVI 

II. DEFINITION OF THE MIXED LATTICE GAS MODEL AND ITS SOLUTION 
ON THE BETHE LATTICE 

We study the grand-canonical version of the mixed lattice gas model defined on a lattice. 
In this model, two types of particles are present, and only excluded volume interactions 
between them are considered, so that the model is athermal. The particles of type 1 (small) 
are such that, when placed on a lattice site, they occupy this single site only, excluding 
other particles from it. Particles of type 2 (large), when placed on a site, occupy its first 
neighbor sites also. Figure [1] shows this model on a square lattice and the generalization to 
other hypercubic lattices is straightforward. While no phase transition is found for the case 
where only small particles are present, it is well established that for a pure system of large 
particles a continuous phase transition in the Ising universality class occurs j^. 

Here we will solve the model with both particles present in the grand-canonical ensemble 
on the Bethe lattice. The parameters of the model are the activities of small particles 
Zi = exp(/ii), where fii is the chemical potential of a small particle divided by ksT, and of 
large particles Z2 = exp(yU2)- We proceed defining the model on a Bethe lattice, which is the 
core of a Cayley tree with general coordination number q. We then consider partial partition 
functions (ppf's) of the model on subtrees with fixed configurations of the root site, which 
may be empty (0), occupied by a small particle (1) or by a large one (2). Considering the 
operation of connecting q — 1 subtrees with a certain number of generations of sites to a new 
root site, we may build a subtree with an additional generation and write down recursion 
relations for the ppf's. If we call gi the partial partition function of a subtree with root site 
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configuration i = 0, 1, 2, we have the following recursion relations for the pff 's: 

g'o = {9o + 9i + 92Y-\ (2a) 
g[ = z,{g, + g,r-\ (2b) 

92 = ^29t^- (2c) 

The prime denotes subtrees with an additional generation. Let us define ratios of the ppf 's 
as Ri = Qi/go, where now the configuration index i assumes the values 1 and 2. From the 
recursion relations for the ppf 's, we may obtain the ones for the ratios, which are: 

R' = Z2-, —-r. (3b) 



The thermodynamic behavior of the model is defined by the values of the ratios after 
a large number of iterations of the recursion relations Eqs. (|3]). We find that, depending 
of the values of the activities, the recursion relations converge either to a fixed point or to 



a limit cycle of period 2. For convenience, we may define the variables x = 1 + i?i and 
1/ = 1 + i?i + i?2. The fixed point equations in these variables will be: 

[x - l)y''-^ - zix'i-^ = 0, (4a) 

{y - x)y'^-' - ^2 = 0, (4b) 



When the recursion relations converge to a limit cycle, the values of the ratios in the core of 
the tree will display a layered structure, so that the ratios in sites in consecutive generations 
assume alternate values. It is convenient, in this case, to define two sublattices (A and B), 
whose sites are associated to the two values of the pair of variables x, y. The equations 
defining this limit cycle values are: 

0, (5a) 
0, (5b) 
0, (5c) 
0. (5d) 



{XA - l)y'ii^ - Zix'^B^ = 

(jjA - XA)y%'^ - Z2 = 

{xb - 1)2/^^ - Zix"^^ = 

{VB - XB)y'A^ - Z2 = 
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Although we were not able to find general solutions for these equations, it is not difficult 
to solve them numerically. This set of equations may be reduced to finding the roots of a 
polynomial in the variable g = x/y, given by: 

h{g) = {l-g){l + z,g^-' + z^r-y-' - Z2f = 0, (6) 

where / = l/y is given by: 

f = g-^(l-g)(l+Z,g^-'y-\ (7) 

The fixed points or limit cycles will correspond to roots of this polynomial in the range 
g G [0, 1] with non- negative values for /. Both stable and unstable fixed points will be 
found. Another numerical procedure to find the thermodynamic properties of the model is 
to iterate the recursion relations Eqs. ([3]) directly, this will lead only to the stable fixed 
points. To study the stability of the fixed points, it is useful to obtain the jacobian of the 
recursion relations. The elements of the 2x2 jacobian matrix calculated at the fixed point, 
ji,jix,y) = dR'JdRj, are: 

(g - l)zix'i-'^{y -x) 
•^1.1 = ' (8a) 

_ {q-l)zix'i~^ 

- 1)^2 s 
32,1 = J2,2 = . (8c) 

The stability limit of the fixed point may then be found requiring the dominant eigenvalue 
of this to have a unitary modulus. The jacobian for the limit cycle J2 will be the product 
of two jacobian matrices defined above, calculated at the pair of variables at the limit cycle, 
so that J2 = 3{xA,yA) X 3{xB,yB)- 

The grand-canonical partition function of the model on the Cayley tree is obtained con- 
sidering the operation of attaching q subtrees to the central site of the tree. If the central 
site is in sublattice A, this leads to the following expression: 

Ya= glBiyB + Zl^B + ^2), (9) 

and a similar expression with the sublattice indexes interchanged is obtained for sublattice 
B. It is easy then to write down expressions for the densities of sites with small and large 
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particles in the center of the tree. They are: 



^1 



Pi,A = ^— g— , (10a) 

P2,A = , (10b) 

Vb + ^l^B + ^2 

(lOc) 

and the densities on sublattice B are obtained permuting the sublattice indexes. 

The free energy in the core of the Cayley tree, which corresponds to the Bethe lattice, 
may be obtained by ageneralization of Gujrati's argument Q, which may be found for a 



particular model in 13|]. Here we present a simple derivation generalizing the one which is 
found in [l^. On the Cayley tree, we admit that the reduced free energy per site (which 
corresponds to the grand-canonical free energy divided by ksT and the number of sites in 
the lattice, which is proportional to its volume) of sites in the m'th generation of the tree 
will be 0*^™'''. For a tree with M generations, starting to count at the surface (m = 0), we 
may then write the total free energy as: 



M-l 

For a tree with one more generation, the free energy in terms of the free energies per site 
will be: 

M 

$(M+i) _ ^(M+i) ^ ^ J2{q - 1)>(^^-^). (12) 

1=0 

Now we notice that: 

_ _ = + q^iM) _^q_ ^)^(M)_ (^3) 

In the thermodynamic limit, we expect the free energies per site to reproduce the sublattice 
structure in the core of the tree, and since sites in consecutive generations belong to different 
sublattices, we notice that in both possible cases (A or B sublattice at the central site) we 
reach the conclusion: 

where the indexes stand for the sublattice. The free energy per site in the core of the tree 
is therefore given by: 

I I y(M+l) 

<^>b = 2 + ^b) = "2 In [y(M)]fa-l) - (15) 
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Using the Eq. ([9]) and the fixed point equations @, after some algebra, we find the result: 



-im 

2 



(16) 



XVA - Xa + VB - Xb + 

We notice that this expression is invariant under permutation of the sublattice indexes. For 
Z2 = 0, where the model is solved trivially, the Bethe lattice calculation furnishes the exact 
solution. The fixed point value in this case is simply Ri = zi, the reduced free energy per 
site will be 06 = — ln(l + zi), and the density of particles becomes: 

zi d(f)b , . 

III. THERMODYNAMIC BEHAVIOR OF THE MODEL 

For simplicity, we will start the study of the thermodynamic properties of the model in 
the limit zi ^ 1, where the density of small particles is very small. In this region, the 
fiuid-solid transition is continuous, and the stability limits of both phases are coincident. 
To obtain the critical line, we may consider the fixed point equations (jl]) and require the 
leading eigenvalue of the jacobian Eq. ([8]) to be equal to -1, so that: 

1 + + J2,2 + jl,li2,2 - jl,2j2,l = 0. (18) 

Now we solve these three equations up to first order in zi, supposing that Z2 = a + bzi, 
X = 1 + biZi, and y = a2 + 62-^1- We are thus lead to the following values of the expansion 
coefficients defined above: a = (q - l)'?"V(g - 2)«, b = -1, bi = [{q - 2)/(g - 
(^2 = {q ~ ~ 2), and 62 = 0. For the particular case zi = 0, this solution has been 

obtained a long time ago [6']. We also notice that the critical line has a negative initial slope. 
This may be understood physically realizing that the presence of few small particles does lead 
to an effective entropic attractive interaction between the large particles, thus favoring their 



ordering [15[. At higher values of zi the slope becomes positive, and finally the transition 
becomes discontinuous, the critical line meets the coexistence line at a tricritical point. 

We may also study the phase diagram in the limit zi ^ 1. In the fiuid phase, we have 
the asymptotic fixed point values x ^ Zi and y ^ Zi. In the solid phase we have xa ~ 1, 
yA ~ 1, Xb ~ ^1, and ys ~ -21 + ^2- From these expressions, we may find the behavior of 
the coexistence line in this limit by requiring the free energies, Eq. (1161) . of both phases to 
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be equal, and this leads us to the asymptotic behavior Z2 ~ z\ for the coexistence line for 
large values of z\. 

The phase diagram in the (zi, Z2) plane is shown in figure [3] for a Bethe lattice with g = 4, 
similar diagrams are found for other values of g > 2. For Zx larger than the tricritical value, 
there is an interval of values of z^ where both fixed points are stable, thus characterizing 
a coexistence of both phases, at a value of Z2 for which both free energies are equal. The 
coexistence line is located between both stability limit lines, as expected. The numerical 
determination of the coexistence line has to be done carefully, particularly close to the 
tricritical point, where the range of values of for which both fixed points are stable 
becomes very narrow. The precise calculation of the localization of the tricritical point 
also demands some care. It is quite easy to calculate the stability limit of the symmetric 
fixed point, since the problem may be reduced to finding the solution of an equation in one 
variable. Once this line is found precisely, we solve the asymmetric fixed point equation on 
it, starting at a value of z\ larger than the tricritical value. As the value of z\ is lowered, 
the largest eigenvalue of the jacobian Ai increases and the values of the variables xa and 
xb which solve the fixed point equations (E]) become closer, as do the variables yA and yB- 
At the fixed point, the dominant eigenvalue is unitary and xa = "^b- This is illustrated in 
figures m thus leading to an estimated position of the tricritical point. Table ([T]) presents the 
locations of the tricritical point for several values of the coordination number g, as well as the 
values of the densities at this point. An alternative procedure for determining the location 
of the tricritical point will be presented below. In all cases we studied, the critical value of 
Z2 for z\ = 0, Z2^c = CL given above, is smaller than Z2 at the tricritical point, therefore for 
Z2 below Z2,c and above the minimum critical value of Z2 the solution displays a reentrant 
behavior as zi is increased, starting in the fiuid phase, then getting into the solid phase and 
finally ending in the fiuid phase again, with two continuous transitions. 

The phase diagram in the plane defined by the densities of small and large particles is 
presented in figure El Again we may obtain the asymptotic behavior in the limits 2:1 ^ 1 
and Zi ^ 1 using the results for these limits presented above. In the limit of low density 
of small particles, we find that the densities are given by pi ~ (g — 2)'^zi/[g(g — 1)'^^^] and 
P2 ~ l/q — {q — 2Yzi/[q{q — 1)'^"^], the critical line for pi ^ 1 shows a linear behavior 
pi ~ 1/g — P2- In the high density limit, we find the following densities on the coexistence 
line: for the fiuid phase pi ~ \ — 1/zi and p2 ~ l/z\~^, so that for this phase we find 
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Pi ^ 1 — p2^^'^~^\ Therefore, the hne reaches the point pi = 1, p2 = with infinite slope. 
For the sohd phase, we get pi = (pi^ + Pib)/2 ^ l/{2zi) and p2 ~ 1/2 — l/(22;i), so that 
Pi = 1/2 — p2. Quahtatively, these features of the phase diagram are similar to the ones in 
Fig. 3 of the paper by Poland jll|. 

In general, the fixed point or limit cycle may be obtained solving equation ([6]) for the 
variable g. For a given value of zi, at small values of Z2 we find a single root for the equation 
in the interval [0, 1]. For values of Zi below the tricritical value, above the critical value of 
Z2, two roots are found, one of them being a double root. They correspond to the fixed point 
variables x and y at sublattices A and B. If zi > zi^tci for ^2 above the limit of stability of 
the solid phase and below the limit of stability of the fiuid phase, five roots are found, such 
that the smaller and the larger ones correspond to the fixed points associated to the solid 
phase, while the intermediate root corresponds to the fiuid phase. Between the extremal 
and the central roots, there are two additional roots, which are unstable. Finally, in the 
region where only the solid phase is stable, three roots are found, the intermediate one being 
unstable. These findings lead to an alternative procedure to locate the tricritical point, since 
there the polynomial h{g), its first and second derivatives vanish. This is discussed in some 
more detail in the appendix. 

Another interesting feature of the Bethe lattice solution of this model is that the isobaric 
curves of the total density of particles as a function of one of the fugacities show a non- 
monotonical behavior in the fiuid phase. Let us define the total density of particles as 
p = pi + 2p2, so that it will be in the interval [0, 1]. If we recall that the grand-canonical 
free energy per site is related to the pressure by: 

^ (19) 



where Vq = V/N is the volume per site, so that we may define 11 = —(pi, to be the reduced 
pressure. For a fixed value of the pressure, we may now obtain the density of particles p as 
a function of the activity zi. Some resulting curves may be seen in figure [61 We notice that 
the curves are not monotonic, starting with a negative slope at 21 = 0. The minima in the 
isobars are located on the dotted lines in the phase diagrams of figures [3] and El We notice 
that this line starts at a particular point of the zi axis and ends at the tricritical point. 
Some aspects of these minima may be discussed analytically. To find the point of the Z2 = 
axis where the minimum of the isobars is located, we may obtain a solution of the model 
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for Z2 <C 1, since we can solve it exactly for Z2 = 0. The fixed point values of the ratios are, 
up to linear terms in Z2: 

R, ^ zi(l- ) , (20a) 



:i + ^1) 

R2 ^ 7 (20b) 

(l + ;zi)'?-i ^ ^ 

Using these approximate solutions and the Eqs. ( ITOj) for the densities, we may then find an 
approximate expression for the total density of particles: 

, o ~ ^1 (g- l)zi -2 
The minimum of the isobars correspond to the condition: 

This derivative may be calculated at ^2 = noticing that: 

dp\ f 9p\ f 9p\ f dz2 



dzijj, \dzij \dz2J ^^\dzij^' ^^^^ 
and since dcpb = —pidzi/zi — P2dz2/z2, the last derivative in the expression above is: 

(^) (24) 

\dzij^ P2Z1 



Finally, we get: 

dp\ {q- l)zi - 1 
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(25) 



so that in the limit Z2 the minimum of the isobars is located at 2:1 = l/(g — 1). A similar 
analysis may be done close to zi = 0, showing that the slope of the isobars is negative there. 
As can be seen in Eq. fl2^ . for constant pressure Z2 is a decreasing function of Zi, a feature 
which is valid in the solid phase also. Therefore, any isobar which starts at {zi = 0,^2), 
will end at {zi,Z2 = 0) such that 11(0, 2:2) = n(-2;i,0). This is apparent in figure |6l where 
we notice that each isobar for finite pressures ends at a finite values of Zi and p. Isobars 
which start in the solid phase, at (0,^2) with Z2 > {q — iy~^/{q — 2^, will cross either the 
critical line or the coexistence line before they end in the fluid phase. In the first case we 
notice a discontinuity in the slope of the isobars. In the second expected the density 

is discontinuous as the coexistence line is crossed, so that the minimum is located on the 
coexistence line. 
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It is worth mentioning that such non-monotonic behavior of the density of particles for 
constant pressure as a function of a fieldhke variable is found in nature, one of the most 
studied examples of this kind is the density anomaly of liquid water, where a maximum is 
found in the isobaric curves of density as a function of the temperature close to the freezing 
pointy^]. In many studies in the literature, simple models were proposed which show such 
anomalies, and it is believed that an interparticle potential with two length scales may be 
the source of the thermodynamic anomalies. A recent work of this kind, where also many 
earlier studies are referenced, may be found in [17]. For such models, the solution on tree-like 



lattices may also be useful IJ], and in this particular example both maxima and minima of 
the isobaric curves of density as a function of temperature were obtained. Although of course 
the present model is very different from the lattice gases related to water, it is interesting 
that here also two length scales are present in the interparticle interactions. 



IV. CONCLUSION 

The Bethe lattice solution of the athermal lattice gas with two kinds of particle, a small 
one occupying a single site and the other a site and its first neighbors, shows a phase 
transition between a fluid phase, at low values of the activity of the large particles Z2 and 
a solid phase, which appears at higher values of and where one of the sublattices is 
preferentially occupied by the large particles. The results of the Bethe approximation for 
the particular case where no small particles are present {z\ = 0) are recovered j^. The 
transition remains continuous for small values of Zi, the critical line starts with a negative 
slope as zi is increased, passes through a minimum and ends at a tricritical point, so that 
for Zi larger than the tricritical value the transition is discontinuous. 

The behavior of the densities of particles, as shown in figure [5l may be compared with 
similar results obtained using series expansions for the model defined on the square lattice. 



shown in figure 3 in reference llj by Poland. Besides the expected quantitative differences, 
we notice in general a qualitative agreement of both diagrams. A significant difference is 
that in our calculation the lines corresponding to the fluid and the solid phases meet at 
the tricritical point forming an angle, while in the diagram by Poland a smooth junction is 
suggested. Since the dotted lines in Poland's diagram are the result of an extrapolation, it 
seems that this aspect may be worth to be studied in more detail using other techniques. 
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However, it may be possible that the angle we find here is a characteristic of the mean field 
approximation implicit in Bethe lattice calculations. 
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Appendix A: Determination of the tricritical point 

As mentioned in the text, one way to determine the tricritical point is to solve the set of 
three equations for the polynomial h{g,zi,Z2) defined in Eq. ([6]): 

0, (Ala) 
0, (Alb) 

0. (Ale) 

Although it seems to be a rather simple task to solve this system of nonlinear algebraic equa- 
tions for g, Zi, and Z2, standard numerical methods, based on Newton- Raphson procedures, 
often do not converge to the expected (physical) solution. This may be due to the fact that 
the first and second derivatives of h with respect to g vanish at the solution. Therefore, we 
used another procedure, taking advantage of the fact that /i is a polynomial in the variable 
9- 

N 

Hg, zi, Z2) = ^ hi{zi, Z2)g\ (A2) 

where N = \ + {q — 1)^[1 + ((? — 1)^]. Now g^, the value of g at the tricritical point, is a 
triple root of the polynomial, so that we may write: 

N N-3 

h{g, zi, Z2) = ^ hi{zi, Z2)g' = {g - gof ^ fiizi, Z2)g' (A3) 

i=0 i=0 

Comparing terms with the same powers of g in the equation above, we may solve for the 
N — 2 coefficients fi in terms of the coefficients hi and gQ. The result is: 

M90, zi,Z2) = - Cj.i (A4) 
13 



h = 
dh _ 
dg 

d^h _ 
dg^ 



for 2 = 0,1,2,..., AT -3, and 



The remaining equations, corresponding to the powers g , g , and (7 , are: 

hN+2 - fN-5 - 3gofN~4 - "iglfN-z = (A6a) 
hN-1 - fN-4 - Sfifo/iv-s = (A6b) 
hN - fN-3 = 0. (A6c) 

The solution of these equations leads to the activities and the value of the variable g at the 
tri critical point. 
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FIG. 1: Two molecules of type 1 (hatched) and two of type 2 (cross-hattched) placed on sites of a 
square lattice. The first may be represented as squares of side a and the second as tilted squares 
of side a\/2, where a is the lattice parameter. 




FIG. 2: Configuration of a model similar that of Fig. [T], which was solved in a particular case on 
the square lattice The possible locations of the large particles are dislocated with respect to 
those of the small particles. 
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FIG. 3: (Color on line) Phase diagram of the model in the plane defined by the activity fractions 
zi/{l + Zl) and 22/(1 + Z2). For a given value of zi, the fiuid phase is stable for lower values 
of Z2, and the solid phase becomes stable at higher values. The transition may be continuous or 
discontinuous, both regimes are separated by a tricritical point, represented by a circle (blue). The 
full line (red) corresponds to the continuous transition and the coexistence of both phases happens 
at the dashed line. In the inset the same diagram is shown with the axes defined by the activities. 
The dotted line (green) corresponds to the minima of the isobaric curves of the total density of 
particles p = pi + 2p2, as discussed in the text. Bethe lattice with q = A. 



q Zl Z2 pi P2 

3 2.74598 11.9231 0.3937(8) 0.1476(5) 

4 1.16956 3.02938 0.2985(7) 0.1117(5) 

5 0.734355 1.539662 0.2407(3) 0.0894(3) 

6 0.533384 0.995756 0.2013(3) 0.0746(3) 

7 0.418259 0.725526 0.1729(2) 0.0641(3) 

8 0.343861 0.566972 0.1515(2) 0.0561(2) 

9 0.291885 0.463718 0.1349(2) 0.0499(2) 

10 0.253531 0.391509 0.1215(2) 0.0450(2) 

TABLE I: Values of the activities and densities at the tricritical point for Bethe lattices with 
different coordination numbers q. The error in the activities values are ±1 in the last decimal 
place. 
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FIG. 4: a) Behavior of the leading eigenvalue of the jacobian of the asymmetric fixed point on 
the limit of stability line of the symmetric fixed point, b) Values of Xa — xb from the solution of 
the asymetric fixed point equation, calculated on the limit of stability line of the symmetric fixed 
point. Bethe lattice with q = 4. 
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FIG. 5: (Color on line) Phase diagram in the plane defined by the densities of particles pi and p2- 
The critical line (red) ends at the tricritical point represented by a circle (blue). The densities of 
the coexisting fluid (black) and solid (green) phases with the same free energies in the two-phase 
region are connected by tielines. Bethe lattice with q = 4. 
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FIG. 6: (Color on line) Total density of particles p = pi + 2p2 as a function of zi for fixed values 
of the pressure 11 (thin full lines) . The starting points of the isobars corresponds to large particles 
only and at the endpoints only small particles are present. The full line (red) corresponds to the 
critical condition and the two dashed lines are the densities of the coexisting phases. For isobars 
which cross the coexistence line, the minimum is located on this line. The loci of the minima of 
the isobars is represented by the dotted line (green). From lower to higher starting densities, the 
isobars correspond to reduced pressures 11 = 0.3, 0.4, 0.5, 0.6, 0.7, 0.86062 (the isobar with the 
minimum at the tricritical point) and 1.0. Bethe lattice with q = 4. 
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